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We develop a finite-temperature perturbation theory for quasi-one-dimensional quantum spin sys- 
tems, in the manner suggested by H.J. Schulz in Phys. Rev. Lett. 77, 2790 (1996) and use this 
formalism to study their dynamical response. The corrections to the random-phase approximation 
formula for the dynamical magnetic susceptibility obtained with this method involve multi-point cor- 
relation functions of the one-dimensional theory on which the random-phase approximation expan- 
sion is built. This "anisotropic" perturbation theory takes the form of a systematic high-temperature 
expansion. This formalism is first applied to the estimation of the Neel temperature of S=l/2 cubic 
lattice Heisenberg antiferromagnets. It is then applied to the compound CS2CUCI4, a frustrated 
S=l/2 antiferromagnet with a Dzyaloshinskii-Moriya anisotropy. Using the next leading order to 
the random-phase approximation, we determine the improved values for the critical temperature and 
incommensurability. Despite the non-universal character of these quantities, the calculated values 
are different by less than a few percent from the experimental values for both compounds. 



I. INTRODUCTION 

Quasi-one-dimensional magnets are notoriously diffi- 
cult to tackle. The backbones of those compounds, 
namely the spin chains, are by now very well understood, 
in some cases even by analytical methods. But until now 
no natural and efficient framework has been developed 
to describe their behaviour when they are coupled by a 
weak interchain exchange J± . 

Useful results have nonetheless been obtained by com- 
bining one-dimensional exact results with a random- 
phase approximation. (RPA) approach to cope with in- 
terchain couplingara. Recently such a method has 
even been applied to frustrated quasi-one-dimensional 
systemsu, yielding sensible predictions. 

From the RPA formalism for the dynamical suscepti- 
bility, one can deduce estimates for nan-universal quan- 
tities, such as the Neel temperatures, or the possible 
incommensurate order developing below the transition 
in a frustrated antiferromagnctu. This is made possible 
by recent progress in the (exact) determination of spin 
chain two-point correlation functions in the low-energy 
regime. The RPA formalism together with those exact 
results are able to cope with exchange anisotropy, and/or 
Dzyaloshinskii-Moriya interaction. 

This approach has been successful in the sense that it 
yields satisfactory results when compared to experimen- 
tal measurement (in some cases even though the inter- 
chain ratio J±/J\\ is not small where J11 is the exchange 
coupling along the easy axis). This owes to the fact that 
on one hand the ratio T c / J» is small enough so that the 
collective one-dimensional excitations have a significant 
influence on the physics at the transition temperature (it 
is always the case for high enough temperature), and on 
the other hand in the case of a cubic lattice, T c /J± is big 
enough. However the RPA is an uncontrolled approxi- 
mation. 

V.Y. Irkhin and A. A. Katanin have calculated correc- 
tions to RPA for spin-1/2 quasi-one-dimensional cubic 



lattices^. Their calculations owe to T. Moriya's empir- 
ical improvement, to the RPA formula for the dynami- 
cal susceptibility^ and it differs notably from what fol- 
lows. Their work has found applications in the estima- 
tion of the Neel temperature of cubic lattice quasi-one- 
dimensional antiferromagnets : KCUF3, S^CuOa and 
Ca2Cu03. Their estimation deviates from the RPA re- 
sult by 25%. 

In the following, we will develop a systematic expan- 
sion and will embed the RPA formula for the dynamical 
susceptibility in it as a natural leading order approxima- 
tion. We will be mainly concerned with lattices made of 
S=l/2 Heisenberg spin chains. YeLthe formalism equally 
applies to anisotropic spin chains! Such an expansion 
has been developed by E. ArrigonU for the physics of Lut- 
tinger liquids. The main differences with our approach 
are the following. Because it is at T = 0, he resums an 
infinite proper set of cumulants. On the contrary, we are 
at finite temperature and we will use the temperature as 
an additional energy scale in the disordered phase. Be- 
cause of this energy scale, we need not resum all of the 
higher cumulants to get a sensible result. But we do have 
to resum temperature-dependent diagrams at the level of 
any cumulant (four-point correlation function for the ex- 
amples given in this work) , because we intend to use our 
expansion down to the critical temperature, where those 
diagrams contribute. 

In section II of this paper, we expose the formal steps 
leading to an extended perturbative expression of the dy- 
namical susceptibility, in terms of a self-energy of the 
two-spin correlation function. 



In section III 



we discuss 

some of the results and peculiarities of this expansion. In 
particular we show that the expansion can be organized 
in terms of the number of RPA-dressed propagators indi- 
rectly related to the small parameter Jj_/T. This prop- 
agator must be regularized and we hint at how it can 
be done. An integral representation of the first correc- 
In section IV, we calculate the correction 



tion is given. 

due to the leading diagram in J±/T and discuss its ef- 
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feet on the physics of cubic lattices and in particular on 
KCuF 3 . In section |v|, we investigate the effect of the 
same correction in a much more involved case, a quasi- 
one-dimensional (or quasi-two dimensional depending on 
the point of view since the interchain coupling is large) 
frustrated antiferromagnet with a Dzyaloshinskii-Moriya 
interaction : CS2CUCI4. 



II. GENERAL PERTURBATION THEORY FOR 
QUASI-ONE-DIMENSIONAL MAGNETS 

We consider the general quasi-one-dimensional mag- 
netic Hamiltonian Ti — Ti\\ + Ti±, where : 



H± = ^Jx^frj)^. 



(1) 



The summation over the spin components is implied 
whereas the latin indices stand for the sites of the 
spins. The quasi-one-dimensional magnetic crystal can 
be viewed as a set of spin chains along which the ex- 
change couplings are supposed to be dominant. H\\ is 
then defined as the part of Ti. which connects spins on the 
same spin chain, whereas Ti± connects by definition spins 
belonging to different spin chains. We aim at giving a sys- 
tematic perturbative expansion of the finite-temperature 
generating functional : 

Z[rf]=Tr T T exp ( — drU- dr^^.S^j , 

(2) 

where (3 — 1/ksT. Now we define the isolated spin chain 
finite-temperature generating functional : 

Z||[ , 0] = Tr T r cxp^— J drTi^—J dr^~] ipj.Sjj 



If we denote Z^ = Z\\ [0], then the average of the observ- 
able O with respect to this functional is 



'MS?]) 



1 

z 



-o 



8tf 



Z 



(4) 



With those notations, we have 



Z[ip] = Zn (exp 



drTi, 



J dr^2 & ■ S^j ] 



(5) 

In a very similar fashion as was done inEl for coupled 
Luttinger liquids, we now introduce a vector field <pi(r) 
in order to perform a Hubbard-Stratonovitch transform 
on Tt± : 



Z[il>] 



Z\\ I T>(f) exp 



1 



(exp (- jf drJ2($i + $i) ■ Sij ) || 



(6) 



The functional integration on (f> corresponds to an inverse 
Laplace transform. The second part of the integrand, 
which is Zu [ip + tft] corresponds to a generating functional 
of the one-dimensional theory with current source ip + 
<t>. Then — In Z\\ [tp + <p] is the free energy of the sum of 

the individual spin chains : In Z\\ [ip + <p] = W[(j> + tp]. 
The summation over the spin chains is included in the 
functional W, which has a Ginzburg-Landau expansion : 

W[$] = ~ f d(l)d(2) Cfl(l, 2) + WS\ > (7) 



where W\ [<p] is the interaction functional : 



1 r 4 



41/11 ••• ' /'.".A.„ ' - - - • : I, r (2)^(3)^(4) 

WH, (8) 



where Jd(i) — f^dri f^° 00 dx^2 n with n the index of the 

spin chain. C' p '(l,..,p) is the time-ordered imaginary- 
time p-point correlation function of an isolated spin 
chain. 

We now work in momentum space and Fourier trans- 
form the functional integrals. We therefore adopt the 
new convention : Jd(i) = ^J^J^^J^^ 
in the case of a three-dimensional magnet. The sum- 
mation indexed by n is performed over the Matsubara 
frequencies u> n — 2im/j3. Finally let us define the field 
theory : 

(0)i = (v$ [exp [F\<t>]) 0\l [ V$ exp (F[<p]) , (9) 



^1)^-1) 



(3) with the weight (<Pm stands for <p a k\)) 



m = - Jd(i) [pjj^w+c^a) 



(10) 



Rewriting Eq. ([]]) in terms of the theory defined by Eq. 
(St), we obtain 



zW\ 
z[6] 



exp 



\ fd(l) [2 Jx]" 1 ,, (1)^)^-1) 



x (exp f- J d(l) [2 Jj-l (1) ^d^-d) {U) 

Interpreting the averaged exponential in ( pd| ) as a 
generating functional and introducing the self-energy 
Y,^ v (w, k) for the two-point correlation function, we de- 
duce that to second order in tp, one has (assuming for 
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simplicity SU(2) invariance so that CjxJ, — C^S^^, 
J^ u = J±_5^ v and T,^ = £<5^„) 

In (Z$\/Z[fy = 0(|^| 4 ) + 

1 + 2Ji(l) (CW(1) + E(1)) ' ^ ( - {)L2) 



5 Ml) 



This form of the two-point correlation function has been 
suggested by H.J. Schulz inH. 



III. CALCULATION OF THE FIRST 
CORRECTIONS TO THE RPA DYNAMICAL 
SUSCEPTIBILITY 

A. RPA formula for the dynamical susceptibility 

To the lowest order of approximation, one can set 
Ti(w,k) = in Eq. jl^). We can then continue analyti- 
cally (on the frequencies) the one-dimensional two-point 
correlation function and therefore recover the dynamical 
magnetic susceptibility 



X3d(w, k) 



Xid(w, k x ) 



1 - 2 J±(k)xid(u,k x ) 



(13) 



As a consequence, the RPA approximation for quantum 
spin systems appears as the leading order of a more gen- 
eral expansion scheme. 



B. Higher-order corrections to the RPA formula 

This clearly shows that a systematic expansion can be 
used. The free (Euclidian) propagator of the effective 
theory is an RPA-dressed propagator (simply called G 
thereafter). Its inverse can be read off from Eq. ( |l0|) : 



[G]- 1 = C& + [2J x y 



(14) 



In real space and imaginary-time, it is given by : 

C(t r\ — — ^ X ^ y ^ Z r ik-r+iw n T 



P^J-oo 2^ Jo 2tt 2tt 

2J±(k) 
X l + 2J ± (fc)C( 2 )(zw„,fc) 



(15) 



Depending on the value of the temperature, the integral 
might be improper and it is then meant that the principal 
value of the integral has to be taken. We postpone the 



discussion on this issue to subsection III D 



The vertices of the perturbation theory are given by 
the multiple 2n-point correlation functions of the spin 
chains. In the case of the spin S=l/2 those are known 
exactly in the asymptotic limit. 

The vertices of the effective field theory (|l0|) involve 
separated space-time points (r, x) and therefore always 



depend on (w, k x ) when written in momentum space. On 
the other hand, they are point-like vertices as far as the 
transverse space coordinates are concerned or, stated dif- 
ferently, do not depend on k± in momentum space. As 
a consequence, all diagrams in the expansion of the self- 
energy T,(u>, k x , k±) are expected to depend on k x . Yet, 
only those with internal RPA-dressed propagator lines 
which are true functions of the input transverse momen- 
tum k± are to depend on it. 

The first fcj_-dependent diagram possesses three RPA- 
dressed propagators and two four-point vertices as de- 
picted on Fig. (|l|). What happens for the transverse 
momenta is reminiscent of what occurs to many-body 
field theories of electron gas, where the dependence on 
the space momenta appears only to the order of this 
diagram, whereas the first diagrams (Hartree-Fock) de- 
pend only on frequencies. For diagrams which do not de- 





(a) 



FIG. 1: Diagram (a) does not depend on k± whereas diagram 
(b) containing RPA internal lines depending on the input k± 
does. 

pend on the transverse momenta, and which are therefore 
merely one-dimensional, one can resort to the simplified 
one-dimensional RPA-dressed propagator : 



G(r, x) = G(t, x, f± = 0) 



(16) 



It is the propagator which has been used inn. 

What is the small parameter of this expansion ? In 
the simplest case of the cubic lattice, it is likely to be 
J±/T. More precisely it is AJ±/T where A is a prefac- 
tor, possibly weakly dependent on the temperature. The 
prefactor A will be given later (Eq. (|l7|)) in the case of 
S=l/2 Heisenberg spin chains. Indeed, each RPA-line 
contributes by an obvious factor of (AJ±/T) 2 in any di- 
agram. Yet in each RPA-line expression remains a non- 
polynomial dependence on AJ±/T corresponding to the 
usual RPA re-summation of transverse paths. Undoing 
this RPA summation, the propagator can be expanded 
in contributions with an exact dependence in (AJ±/T) 2 , 
(AJ±/T) 3 , etc. For a diagram with p RPA-lines in it, it 
is rather (AJ ± /T) 2 p, {AJ ± /T) 2 p +1 , etc. 

So whatever the subtle dependence of the RPA-dressed 
propagator on the temperature, this expansion can gen- 
uinely be seen as a high-temperature expansion in the 
parameter (AJ±/T) 2 . More formally, it is also an ex- 
pansion in the number of RPA-dressed propagator lines 
although their dependence in the small parameter is more 
intricate. 

As a consequence, the conditions of applicability of 
this perturbation theory are that (AJ±/T) 2 <C 1 but 



also T I Jii <C 1 in order for the field-theoretic tools to be 
valid (in particular in the calculations of the spin correla- 
tion functions at finite temperature). For the compounds 
studied here those conditions (which have to be modified 
in the case of a frustrated magnet) turn out to be satis- 
fied. 

The expansion also depends on the dimensionality of 
the lattice. This dependence is obvious at the order of 
RPA, where the small parameter is there .proportional to 
the transverse coordination number (seelj). The depen- 
dence is far less clear at higher order, where the dimen- 
sionality is encrypted in multidimensional integrals. It is 
nevertheless possible to take the d — > oo limit in these 
integrals in order to study this dependence. But this is 
beyond the scope of this work. 



C. Details for the first correction 

Let us take into account the very first correction to 
the dynamical susceptibility. So we consider the first 
non-trivial term in the perturbative expansion of the self- 
energy. It involves the four-point correlation functions 
of the spin S=l/2 Heisenberg spin chain. We decide to 
truncate the Landau-Ginzburg expansion of W\ [4>] to the 
quartic term in ip (six- and higher-point correlation func- 
tion do not contribute at this order anyway). The field 
theory expansion formally resembles a 3-component </> 4 
theory. In particular the very first correction to the self- 
energy is given by Hartree-Fock diagrams (Fig. (|)). The 
"free" propagator of this </> 4 field theory Eq. (|15) is built 
on the usual imaginary-time two-spin correlation func- 
tion of the Heisenberg chain but dressed by the RPA 
corrections. It is therefore a significantly enhanced prop- 
agator and the first correction to RPA in this scheme is 
expected not to be negligible. 



i,.,.,.: 



FIG. 2: Hartee-Fock diagrams with symmetry factors which 
are the first non-trivial terms of the self-energy S(oj, k) 

The only vertices of this truncated field theory are 
given by the four-point spin correlation functions. The 
staggered four-point correlation functions at the isotropic 
point can be computed thanks to bosonisation. We 
will denote by A(A/T) the product of the Lukyanov- 
Zamolodchikov prefactor with the logarithmic correc- 
tion induced by the marginally irrelevant current-current 



correctio 



.4 



A 



A 



(17) 



For clarity, we decompose the imaginary-time four-point 
correlation function into A(A/T) and the purely confor- 
mal part of the correlation functions C X xxx and C X xyy : 

T r( S (0) S (zi) S (z2) S (z 3 )) = A 2 (j^J C xxxx( z l; z 2,Z 3 ), 

rp I QX QX qV QV \ 

J t^ ( o)<5( Zi )J (z2) <3 (z3 )/ 



A 



where C XX xx as well as C XX yy are given by : 



|e(*i)e(ar 2 -z 3 )| 



\e( Z2 )e(z 3 )e{z 1 ~ z 2 )o( Zl - z 3 )\ \e{ Zl )o(z 2 - z 3 )\ 

|e(z 2 )6(zi-z 3 )| 2 



16(^)903)6(0! - z 2 )0(z 2 - z 3 )\ |e(z 2 )e(z! - z 3 )\ 
\Q{z 3 )Q{ Zl - z 2 )\ 

' \e{z 1 )e{z 2 )e(z 1 - z 3 )e(z 2 ~ z 3 )\ 

2 



|9(z 3 )e(zi-z 2 )| 



(19) 



and 



C^ YY (z 1 ,z 2 ,z 3 ) = (-ir +X2+X3 - 



9(zi)e(jS 2 - Z 3 )| 



x Re 



/e(z 2 )9(z 3 )9(zi-z 2 )9(zi-Z3) 
9(z 2 )9(z 3 )9(zi - z 2 )Q(z 1 - z 3 ) 



1 (20) 



where we have denoted 

ii / j^ r J 1 \ 

Q(z = x + ir) = sinh — (x + iur) . (21) 

ttT \ u J 

Gnly the staggered part of the correlation functions 
which dominate has been taken into account. 

The non-universal constant A is taken to be A = 
24.27J|| as calculated irfl. Although the value of A 
and of A are ^somehow different from those extracted 
from numericstSI and used inQ, there is no contradiction 
with the numerical estimates of the correlation functions 
themselvesH. In particular, the Neel temperature esti- 
mated through RPA for KCuF 3 is very close. Discrep- 
ancies might nevertheless appear for a different range 
of temperature and when the self-energy corrections are 
taken into account. Finally u = \J\\ is the spin-1/2 
Heisenberg chain velocity given by Bethe Ansatz. 

The first contribution to the self-energy is then : 



5 



1 /A\ f 13 f 00 
(u>, k) — -A 2 I — J J dTid,T2d,T3 J dx\dx2dx 3 e 



G(z 3 - z 2 ) 



1(7(4) . 



1^(4) 
2 ^xxyy 



(zi,z 2 ,z 3 ), (22) 



where the integrals are performed in real space like was 
done inu. We are mainly interested in the knowledge of 
Y>(uj,k) around k x — tt because the isotropic correlation 
functions are most significant at this point. More pre- 
cisely, most of the spectral weight remains at this point. 
It turns out that, for our purpose, the numerics are much 
more efficient when done in momentum space : 

•A\^ r dq x 



£« k,fc 



1 

W 



T 



E 



2tt 



C {A) 



i c (4) 
2 w xxyy 



(k, k, q n , q„ 



2tt 



dq± 



■G(q n ,gi(fe3) 



where k = (w, k x ) and q„ = (cj n , q x ). 

Now only C^(uj n ,k x ), which enters in G(q„,(f^) 
through Eq. (|l5|), remains to be known. It is the 
spin-spin time-ordered imaginary-time (isotropic) corre- 
lation function. The large distance behaviour of the 
finite-temperature correlation function can be deter- 
mined by combining results obtained from the Bethe 

t-SL^l/2 spin chain, 



Ansatz solutioifflL3 of the Hei 
with field theoretic techniques' 1 

A(A/T) 



ttT/i 



I sinh 



-(x + iut)) 



(24) 

This result can be extended to the anisotropic spin chain 
as well, although A is known exactly only at the isotropic 
point. The frequency and momentum dependence is ob- 
tained by Fourier transformation and analytic continua- 
tion of the time-ordered jjaagjaary-time staggered corre- 
lation function @ (sccfiySEll) : 

j — uk x \ 



Xid(w,7r + k x ) 



A(A/T)T(l 



4ttT 



4T 



r(i- 

J 4ttT 



■ u) — uk x 
' AttT 



r( 



• uj+uk x 
' 4ttT 



(25) 



But we will mostly use this result in its Euclidian form 
(before analytic continuation) : 



-Xid(iu n , k x ) ■ 



(26) 



Finally let us mention the fact that as claimed before, 
£W(k, k±) does not actually depend on k±. We also 
note from the above result that the obvious prefactor of 
JlEW is expected to be (AJ ± /T) 3 . 



D. Prescription for the RPA-dressed propagator 

The RPA-dressed propagator in momentum space (|l4| ) 
may, for some values of the variables (u) n , k) and param- 



eter (T), exhibits a singularity. In the explicit formula 
for the diagrams, this propagator is always integrated 
over and the principal values of the resulting integrals 
are finite. Still, the presence of this singularity has to be 
understood and a correct prescription for it (here taking 
the principal value of improper integrals) to be justified. 
The difficulty arising from its presence can be overcome 
as follows. 

Let us assume that this field theory has a critical tem- 
perature T c , which is the exact theoretical estimate of the 
Neel temperature. The perturbation theory is expected 
to be valid for high enough temperatures and only for 
temperatures above T c in the disordered phase. The RPA 
approach provides an estimated critical temperature T' c 
presumably larger (as will be verified later) than T c . It 
corresponds to a pole in the dynamical susceptibility. For 
T > T' c no singularity is expected to appear in the RPA 
propagator and no problem occurs in the perturbation 
expansion. Whereas when T c < T < T' c , the denomina- 
tor of the propagator is negative which is tantamount to 
realize that subdiagrams given by an RPA line just add 
up to infinity. Nonetheless this is an unphysical singu- 
larity which can be cured by a proper prescription as we 
are going to hint at. 

Let us show how it is done in the case of the calcula- 
tion of T", the critical temperature at the next leading 
order like was done above when calculating S*- 1 ) . We may 
expect T" to satisfy T c " < T' c and hence yield a prob- 
lem. This statement is based on the fact that including 
the effects of the four-spin correlation functions (in ad- 
dition to the two-spin correlation functions) amounts to 
take into account quantum fluctuations more precisely, as 
compared to merely restricting to the Gaussian fluctua- 
tions of the RPA approximation. As we have seen 
depends on a singular RPA-dressed propagator. Now, 
we will add to (which sum will hence be denoted 
SW) subdiagrams which would normally be appearing 
at higher order in the expansion. The single RPA line 
drawn in the diagram for £W appears now as a skeleton 
line in the diagram for S^ 1 ). In this case, it stands for an 
RPA line plus the self-energy correction itself. Y,^ 
is therefore defined as the sum of all irreducible bubble di- 
agrams as drawn in Fig. (||) . It satisfies a self-consistency 
Dyson equation given by the expression ( |23| ) but where 
G(q ni q±) is now replaced by G(q„,<fj_) defined by an 
enhanced version of Eq. (Q) 



(27) 



By construction, the singularity is now avoided, the pole 
of G being displaced to a boundary of the integration 
domain. At this boundary the singularity is genuinely in- 
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FIG. 3: Bubbles diagrams that add up to form S' 1 - 1 



tegrable provided the lattice is at least three-dimensional 



(which is fully consistent with the Mermin- Wagner the- 
orem). Hence subdiagrams have ultimately canceled the 
singularity encountered above. Yet our method is a sys- 
tematic one and does not rely on a self-consistent ap- 
proach. But we can now extract from this construction 
the part of EW which corresponds to the regularized £W 
we have to calculate in our scheme. In order to do so, we 
decompose the integral J * rJl)2 G(q n ,gl) appearing in 
the expression of 



J 



L (^ G(q "^ } 



2tt 



U (27T) 

2 \dq ± 

o 



(2tt) 2 



2Ji 



2Jx (C( 2 ) 

2J ± 
f 2 J ± C( 2 ) 

2J± 



(2tt) 2 1 + 2 J ± C<? 



£« 



(2J ± ) 2 £« 



(1 + 2 J ± (a( 2 ) 



£«)) (1 + 2J X C(2)) 
(2Ji) 2 E« 



(2^) 2 (1 + 2J ± (C< 2 ) + £(!))) (1 + 2J ± C( 2 



(28) 



The first term of the last r.h.s. is the regularized ex- 
pression for the RPA-dressed propagator we use in 
whereas the second term (which also ought to be regu- 
larized) is a correction to it at a higher order in A J±/T. 
The prescription scheme consists therefore (at least in 
this case) in taking the principal value of the integral 
(symbolized by P) which turns out to be finite. 



IV. APPLICATION TO CUBIC S=l/2 
ANTIFERROMAGNETS 

Corrections to RPA in the framework of cubic S=l/2 
antifcrromagnets was the subject ofcl. In particular, the 
authors have derived the integral expression for the di- 
agrams of Fig. (J2J) . However their correction does not 
correspond to the first self-energy correction but rather 
to a subset of diagrams. The self-energy correction de- 
rived here includes the re-summation of the one-particle 
reducible diagrams made of chains of their contribution. 

Let us apply our formalism, using the E' 1 ' correction 
to the self-energy to the compound KCuF 3 . The experi- 
mental value is T c = 39 kM The RPA Neel temperature 
is estimated to be T c — 52.3 K. Taking into account their 
correction, V.Y. Irkhyn and A. A. Katanin then deduced 
T c = 36.7 K, which correction is of order 30%. The singu- 
larity of the RPA-dressed propagator is removed by, using 
a semi-empirical approximation due to T. Moriyatl. But 
its non-trivial dependence on the temperature and the 
couplings constant disappear as well. 

From their calculation, one can deduce the value of 
some intermediate integrals to be calculated. Making use 
of the values of these integrals, and therefore resorting to 



Moriya's approximation, but within the self-energy cor- 
rection approach on obtains the value T c — 31.2 K (that 
is, at this order, re-summing the reducible diagrams). 
This correction of order 40% is significantly stronger. 

To determine the critical temperature with the method 
developed here is somewhat more complicated. Indeed, 
the self-energy is also T-dependent in a non-simple way. 
We have therefore to solve the problem by iteration on 
the value of the estimated critical temperature. The for- 
mal calculations detailed above can be applied with the 
transverse lattice structure factor 



(29) 



J±(k±) = J±(cos(k y ) +cos(fc z )) 



of three-dimensional cubic lattices. From Eq. (|1^ 
and after analytic continuation, we obtain the three 
dimensional dynamical magnetic susceptibility 



X3d(k,fc±) 



xS(k) + £ xx (k,£u 



l-2J x (fcx) Xw( k ) + Sxx ( k . k 



(30) 

The instability condition which can only be satisfied at 
zero frequency is therefore 

2 Jx(£)X(0,fc) = 1, (31) 

where X(k) = XidO 4 ) + Because at this order 

S xx (k, k±) does not depend on k± and because it does 
not a priori change the monotony of x xx (k) with respect 
to k x , one can first maximize the l.h.s. on k. It leads to 
Neel order in the three directions k x = k y = k z — tt so 
that the instability condition is reduced to 4 J± X(0, 0) + 
1 = 0, where J± — J y = J z . 
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For KCUF3, the exchange values are J\l^ 406 K and 
J± = 19.1 K (5 % of the main coupling)o. The small 
parameter close to the transition is A J±/T ~ 0.3. The 
numerical result of this calculation is T c = 40.3 K, fairly 
close to the experimental value. Finally let us mention 
that not taking into account the log-log correction would 
have led us to T c — 37.7 K. So the subtle log-log cor- 
rection would presumably be more significant than the 
second order correction E' 2 '(k, k±). 

On Fig.(^J), we have drawn the general curve of the 
estimated critical temperature T c of cubic lattices as a 
function of the interchain exchange J = J± = J y = J z . 
The upper curve (RPA) corresponds to the RPA estima- 
tion of the critical temperature. The lower curve (IK) 
is deduced from Irkhin and Kataniu,' estimation. It can 
be deduced from their main resuhu reformulated with 
our notations and from the use of the exact correlation 
function prefactor 



T c = k J±A 





rr(i/4)i 


(1) 


|r(3/4)J 



1 2 



(32) 



where k ~ 0.70. The intermediate curve (NLO) corre- 
sponds to our next-leading-order estimation of the crit- 
ical temperature. It is significantly lower than the RPA 
one as expected (about 25%). 




FIG. 4: RPA, next leading order and Irkhin and Katanin' 
estimations of the critical temperature T c of cubic lattices in 
units of Jy as a function of the interchain exchange J = Jj_ = 
Jv = Jz in units of J« . 



an ordered phase. The order is cycloidal. Its parame- 
ter is incommensurate and measured to be ko = 0.186. 
Remarkably, in this phase as well as in the, disordered 
phase, the excitations spectrum is incoherentEllEIBlEi 

Its magnetic Hamiltonian has been recently experi- 
mentally determined with great accuracyEj. It can be 
decomposed as : 



U 



n 



(k) 

plane 



(A;) 
plane 



n 



(k,k+l) 
intcrplane 



J 



(fc,fc+i) 



intcrplane 



E 
E 
E 



s, 



D- S, 



X S; 



i+l,j,k 



(33) 



In particular a Dzyaloshinskii-Moriya interaction (7i DM ) 
has been proven to exist on the interchain exchange 
paths, revealed by its anisotropic nature. Experimen- 
tal estimates for the exchange couplings in CS2CUCI4 are 
J|| = 4.34 K, J± = 1.48 K (about one third of the main 
coupling), J z = 0.20 K aud^nally \D\ = 0.23 K (about 
5% of the main coupling)EJ'Efl. 

So it appears that this compound is essentially two- 
dimensional. One of the two-dimensional spin lattices is 
represented on Fig. (|5|). Although the interchain cou- 



(1) 



J// 



(2) 



(1) 



FIG. 5: Exchange paths within the planes : solid lines denote 
the strong exchange Ju, dashed lines the weaker, frustrated 
exchange J±. Dzyaloshinskii-Moriya exchange of magnitude 
D also stands on the interchain paths. Two types of sites 
among the four in an elementary cell of the lattice are distin- 
guished. 



V. APPLICATION TO FRUSTRATED S=l/2 
ANTIFERROMAGNETS 

A. RPA approach for the dynamical susceptibility 
in the disordered phase of CS2CUCI4 

CS2CUCI4 is a spin- 1/2 frustrated antiferromagnet. At 
a temperature of T c = 0.62 K, it shows a transition to 



pling J±/J\\ is considerable, the smallness of the ratio of 
transition temperature to bandwidth T c /t:J\\ w 0.05 in- 
dicates that the quasi-one-dimensional approach we are 
advocating might be tried on CS2CUCI4. 

An RPA approach has been proven to reproduce qual- 
itative features of the compound (incommensurability, 
asymmetry of the dynamical structure factor around 
k x = ff. as well as by construction, the incoherent 
spectra)H. It also gives very reasonable estimates for the 
critical temperature, and the incommensurability. 
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Given the latest experimental data, and also taking 
into account the subleading correction in the logarithmic 
dependence of the spin correlation functions (which was 
intentionally neglected irn), RPA gives fco = 0.197 and 
T c = 0.85 K. 

Those numbers are off by less than 25% compared to 
the experimental values, which is fairly good given the 
complexity of the physics displayed by CS2CUCI4. It was 
shown that some qualitative features of CS2CUCI4 were 
without doubt reproduced by the approach advocated 
irn. Knowing the corrections to the RPA result would 
hopefully decide whether such a perturbative approach 
can be used on this compound for determining quantita- 
tive results. 

We now briefly reproduce the calculation leading to 
an RPA-formula for the dynamical susceptibility as the 
calculation of the subleading order makes use of it. 

The elementary cell of the CS2CUCI4 crystal has four 
Copper ions. Besides, the vector of the Dzyaloshinskii- 
Moriya (DM) interaction which lies on the interchain 
bounds, points in the third direction (denoted Oz here). 
The orientation of the DM vector is staggered from one 
plane to another : D = ±D e z . 

Hence, the magnetic Hamiltonian is now anisotropic. 
Those complications make the RPA-dressed propagator 
possess a matrix-form which can ultimately (in the com- 
plex representation (H — ) of the quantum spins) be re- 
duced to a 4 by 4 matrix. Eq. ( |l4| ) is still valid under its 
matrix form : 



[G" 1 



: (2) 



[Jj 



(34) 



The four components of the vectorial space in which it 
is defined correspond to the — component of the four 
distinct spins in an elementary cell on the right (resp. + 
component of the spins on the left). 

I4 and 



We have C^l = C^l 



where 



J+K I 
J+K I 

1 J-K 
I J-K 



J(fc) = J± (cos(ky) + cos(k x — k y )) 



(35) 



K(fc) = D (sin(fcy) + sin(k x 
\(k) = J 2 cos(fc z ). 



(36) 



The transverse RPA-dressed propagator of the effective 
field theory is related to the transverse RPA imaginary- 
time correlation function G+^(fc x , u>) (itself related to 
the transverse RPA dynamical susceptibility xtdi^x, 
by analytic continuation) by 



G+- 



+ Jt"G+- =I 4 . 



contributions from the various sub-lattice correlators, i.e. 
by taking e.g. 



£(s+M)s ( - (-^,-fc)>, 



(38) 



where the summation is over the four types of sites. After 
analytic continuation on the frequencies, we obtain the 
following RPA expression for the transverse dynamical 
susceptibility 



xtd ( k > fc ±) 

X+-(k)(l + N 1 (fc) X + -(k) 



l-2J(fc)x+-(k) + N 2 (fc) [x+-(k)]' 
where k = (u>, k x ) and 

Nr(fe) = 1(A)- J(fc) , 

Na(fc) = J 2 (fc) - K 2 (fc) - l 2 (fc). 



(39) 



(40) 



The RPA critical temperature as well as the incommensu- 
rability are then determined through the instability con- 
dition obtained by annihilating the denominator of Eq. 

(j(fc) ± yjwfi) + |2(fc)) X + ~(0,M = 1. (41) 

The relevant instability corresponds to the higher possi- 
ble temperature. In order to solve Eq. (|l]) for it, one 
can maximize the l.h.s. of Eq. @ over k. Then one 
deduces that the instability occurs along the chains di- 
rection k y — k x /2 and that k z = (Nccl order in the 
third direction). Then k x and T c have to be determined 
numerically. 

We have assumed that the main instability is given 
by transverse excitations. So we need not calculate the 
longitudinal RPA propagator to calculate the RPA in- 
stability condition. However, it participates to the next- 
leading-order correction and we shall need it later. Eq. 
( |l4| ) is still valid under its matrix form : 



,(2) 



[2. 



(42) 



The four components of the vectorial space in which it is 
defined correspond to the z component of the four dis- 



tinct spins in an elementary cell. We have 
and 



J I 
J I 

1 J 
I J 



— ^ z 



(43) 



(37) 



B. First correction to RPA 



The transverse time-ordered imaginary-time two-point 
correlation function of spins is obtained by adding the 



In the complex spin representation which is more con- 
venient in the case of CS2CUCI4, the (matrix) dynamical 
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magnetic susceptibilities are after Eq. (11) : 
X+T(k) + E+-(k,fei 



*3d ( k > fc jJ 

X&Mx) 



i-j+-(fc)(x+-(k) + E+-(k,fc L )) 

X f d (k)+E zz (k,fc ± ) 



1 - 2 



JI Z (£) (x?t( k ) + e»M±) 



(44) 



Because of the staggering of the DM vector from one 
plane to another, there is no chirality on the three- 



dimensional spin correlation functions and xtd ( k ' ^x ) = 
X3 d + (k,fc x ). 

In the particular case of the first correction (Hartree- 
Fock-like correction), the matrix EijT appears to be di- 
agonal. Indeed, it is made up of a single four-point cor- 
relation function which involves four spins belonging to 
one type of sites (among four). This does not hold at 
higher order. As a consequence, we may see E^7 as a 
number which is given by : 



I ki i 



A 2 

ft 



E 



dq x 
. 2ir 



C (4) 



Cxxvy ) (k, -k,q„,-q„) 



dq± 
„ (2tt) 2 



[2G+-] (q n ,5l) 



" ft 



(45) 



The three terms in (k, k±) are derived from the 
diagrams on Fig. (^3j). The integrals over k y and k z 
are performed over an extended Brillouin zone (from 
(ky,k y ) € [0,7r] 2 to (ky,k y ) £ [0,2tt] 2 ) and the propa- 
gators expressions below take this extended scheme into 
account. A similar expression can be obtained for E zz ^ 
but is useless for our purpose. 



+ zv 



FIG. 6: Hartee-Fock diagrams with symmetry factors which 
are the first non-trivial terms of the self-energy E" 1 (ui,k). 
The full line propagator is associated with G + ~ , whereas the 
dashed-dotted line corresponds to G zz . The full line box is as- 
sociated with CxtL 
responds to Cxxyy 



whereas the dashed-dotted line box cor- 



The RPA-dressed propagators G + (q n ,sj_) 
G zz (q n , q±) can be derived from Eq. (|34| ) and Eq. 



and 



G+- = 



J + (J 2 



K 2 



l 2 )G+ 



(2) 



2JC ( + } - 
2 (J + l) 



(J 2 - K 2 



l 2 )[G,+ 



(2) 



1 + 2 (J + I) C, 



(46) 



(2) 



where C z ^ = Cf^. /2 is the time-ordered imaginary-time 
spin-spin correlation function of the isolated Heisenberg 
chains. When the temperature approaches the theoreti- 
cal critical temperature, we expect the second contribu- 
tion to the self-energy correction E^y (which depends 



on G zz ) to be quantitatively much smaller than the first 
contribution depending on G H . Indeed the RPA propa- 
gator G zz has an RPA critical temperature (its pole in T) 
much higher than the one for G^ . It is therefore much 
less singular that the latter in the temperature range of 
interest. This has been checked numerically. 

In the case of such a frustrated system, it is less clear 
what the small expansion parameter is. However, we can 
get a rough idea on inspecting the next-leading-order cor- 
rection. First of all there is a prefactor (A/T) 2 , given by 
each RPA-line. In addition there is a dimensionful contri- 
bution coming from the transverse lattice structure factor 
which depends on the interchain exchange couplings. But 
contrary to cubic lattices, it cannot be meaningfully ex- 
tracted from the integral. The RPA-dressed propagator 
G" 1 appearing in the integral is more singular at (T c , ko) 
estimated thanks to RPA. The integral value will there- 
fore be dominated by the value of the integrands when 
k ~ fco, when T is close to T c . So, at least when T 
is close to T c , the expansion parameter is of the order 
of (A J±(ko)/T) 2 . As a consequence, in the case of frus- 
trated quasi-one-dimensional magnet leading at the tran- 
sition to an incommensurate order, A J±/T does not nec- 
essarily have to be small provided A J±(ko)/T is. Note 
that the small parameter in the RPA formula is genuinely 
AJ ± (k )/T. 

After analytic continuation, we obtain the following 
expression for the transverse dynamical susceptibility : 



Xtd (k,fc±) 



(k) (l + N x (fc)X(k) 



l-2J(fc)X(k) + N 2 (fc) [X(k)f 



(47) 

withX(k) = xtd (k)+^m (k). The instability condition 
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at this order is therefore 

J(fc) ± ^K 2 (fc) + I 2 (fc)) X(0, k x ) = l. (48) 

Because at this order S(k, k±) does not depend on k±, 
it is as easy as in the RPA case to maximize the l.h.s. 
on k±. It leads again to a longitudinal instability and a 
Neel order from one plane to another. 

C. Numerical results 

The numerical computations performed to evaluate the 
critical temperature T c as well as the incommensurability, 
are more involved than in the cubic lattice case, where it 
is obvious that a Neel order prevails below T c . The insta- 
bility condition has to be solved with respect to T c and 
fco. The self-energy correction itself, once the obvious di- 
mensionful prefactor has been put aside, depends on the 
temperature T through the RPA propagator and depends 
on the ratio fco /T through the four-point correlation func- 
tion. An iterative algorithm on (T c , kg) can nevertheless 
be used. The small parameter of the expansion close to 
the transition is (-J + VK 2 + I 2 ) (tt + fc ) x A/T ~ 0.4 . 

Our findings are the following. The critical tempera- 
ture is estimated to be T c = 0.66 K to be compared to 
the experimental result T c — 0.62 K. The incommensura- 
bility is estimated to be ko = 0.182 to be compared to the 
experimental result fco = 0.186. This is quite remarkable 
since the errors for the results obtained are less than a 
few percent. 

This in return validates the rougher estimates from 
RPAEl which were already quite satisfying. It makes it 
improbable for the success of RPA applied to quasi-one- 
dimensional magnets to be merely due to chance. 

VI. SUMMARY AND CONCLUSIONS 

We have shown that recent one-dimensional exact re- 
sults from integrable models and quantum field theory 



can be applied to one-dimensional spin- 1/2 antiferromag- 
nets to compute quantities such as critical Neel temper- 
atures. Their computation can be made systematic in 
perturbation theory. On rough grounds, it can be seen 
as a high-temperature expansion in Jj_/T. To the next- 
leader order, the leading order being the random-phase 
approximation, the errors committed differ by less than 
10% from the experimental values at least in the two 
cases investigated above. Although those observables are 
non-universal, the calculation only depends on the mag- 
netic Hamiltonian, i.e. exclusively on the knowledge of 
the exchange couplings. Even incommensurate order pa- 
rameter can be accurately determined this way. At least 
this has be shown on the case of the frustrated compound 
Cs 2 CuCl 4 . 

The perturbation theory allows more generally to give 
a perturbative estimation of the three-dimensional dy- 
namical susceptibility. But it could as well be used to cal- 
culate corrections to multi-spin three-dimensional corre- 
lation functions starting from the one-dimensional func- 
tions. 



Yet it is somehow hazardous to go beyond the next- 
leading order mainly used in this work. Although the 
numerical calculations for higher-order correction are 
achievable, the resulting correction is likely to be within 
the field theory approximation error range. The spin 
two-point correlation functions are indeed only asymp- 
totically exact. Being more precise would require to go 
beyond the knowledge of the (mathematical) equivalent 
of the correlations at large distances. For example, one 
could include-the space-dependent renormalization group 
correctionsEjO. 
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